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Abstract 

It is known that a smooth LTB model cannot have a positive apparent central acceleration. Using 
a local Taylor expansion method we study the low-redshift conditions to obtain an apparent negative 
deceleration parameter q app (z) derived from the luminosity distance D^{z) for a central observer 
in a LTB space, confirming that central smoothness implies a positive central deceleration. Since 
observational data is only available at redshift greater than zero we find the critical values of the 
parameters defining a centrally smooth LTB model which give a positive apparent acceleration at z > 0, 
providing a graphical representation of the conditions in the q$ , q^ pp plane, which are respectively 
the zero and first order terms of the central Taylor expansion of q app (z). We finally derive a coordinate 
independent expression for the apparent deceleration parameter based on the expansion of the relevant 
functions in red-shift rather than in the radial coordinate. 

We calculate q app (z) with two different methods to solve the null geodesic equations, one based on 
a local central expansion of the solution in terms of cosmic time and the other one using the exact 
analytical solution in terms of generalized conformal time. 
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I. INTRODUCTION 



High redshift luminosity distance measurements require a dominant dark energy component 
if interpreted under the assumption of inhomogeneity and isotropy of FLRW cosmological 
models. Because of the unknown physical nature of dark energy and the problems in the 
disagreement between the quantum field theory prediction of the value of the vacuum energy 
compared to the value of the cosmological constant inferred from these observations, some 

nn 

alternatives have been proposed, such as the possibility that we may be located [H, |2[ at the 
center of an inhomogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) 
solution of Einstein's field equations. A general approach to map the luminosity distance as a 
function of redshift Di(z) to LTB models has been developed in 3], showing that an inversion 
method can be applied successfully to reproduce the observed D L {z). Interesting analyses of 
inhomogeneous cosmological models and other related problems are given for example in jj- 



34j. The idea to use galaxy number counts [35] to distinguish between inhomogeneous models 



and ACDM has been investigated both analytically and numerically 36H38] , showing how not 
centrally smooth LTB models can mimick the effects of a cosmological constant both for the 
redshift spherical shell mass mn(z) and Dl(z). More recently there has been clear evidence 
that LTB solutions cannot provide a fully consistent cosmological model compatible with all 



available observations 
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40] . They nevertheless provide a good toy model to understand the 



effects of inhomogeneities which motivate the present study. 

In this paper we compute the low redshift expansion of the apparent deceleration parameter 
q app (z) for an observer located at the center of a centrally smooth LTB model. After confirming 
that the central value of the apparent acceleration q app (0) is positive, we derive the general 
conditions for a negative q app (z n ) at some arbitrarily low red-shift z n . This is motivated by the 
fact that observational data is only available at redshift greater than zero, but still sufficiently 
small to make the Taylor expansion accurate. 

We give a graphical representation of the conditions for a negative apparent deceleration in 
the (<7o PP , Qi PP ) plane, which are respectively the zero and first order terms of the central Taylor 
expansion of q app (z). We calculate q app (z) with two different methods to solve the null geodesic 
equations: one based on a local central expansion of the solution of the Einstein's equations 
in terms of cosmic time, reported in the Appendix, and another one using the exact analytical 
solution, confirming the two methods are in agreement. 

The expansion of the solution in terms of cosmic time is quite useful also for other applica- 
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tions requiring a foliation of space-time in space-like hyper-surfaces, such as spatial averaging, 
which is much more difficult to study using the analytical solution in terms of generalized 
conformal time coordinate. 



II. THE LEMAITRE-TOLMAN-BONDI (LTB) METRIC 



Lemaitre-Tolman-Bondi solution can be written as 4lM43] 



ds 2 



-dt 2 + 



(R, r f dr 2 



+ R 2 dQ 2 



1) 



1 + 2E 

where R is a function of the time coordinate t and the radial coordinate r, R = R(t,r), E is 
an arbitrary function of r, E = E(r) and R, r = dR/dr. 
The first integral of Einstein's equations give 

2 2E(r) + 2M(r) ^ 



R 

R. 



p{t,r) 



R 2 

2M, r 

R 2 R, r 



R 3 



(3) 



with M = M(r) being an arbitrary function of r . We are denoting with a dot the partial 
derivative with respect to t, R = dR(t,r)/dt. There analytical solutions of Eq. fl2]) expressed 
in terms of a time variable r = J** dt'/R(t', r) is 

M(r) 



Y(r,r) 



tT,r) 



-2E(r) 

M(r) 
-2E(r) 



I - cos ( ^-2E(r)r 
1 

sin 



-2E(r) 



-2E(r)r 



+ t b {r) 



(4) 
(5) 



where Y has been introduced to make clear the distinction between the two functions R(t, r) 
and Y(r,r) which are trivially related by 



R(t(r,r)) = Y(r,r) 



(6) 



and tfe(r) is another arbitrary function of r, which plays the role of a functional constant of 
integration and is commonly called the bang function, corresponding to the the possibility of 
inhomogenous big-bang/crunches times. 
After introducing the functions 



A(f,r) = ^2 k(r) 



2E(r) 



Po(r) 



6M(r) 



(7) 



so that Eq. (CQ) and the Einstein equations (12]) and (|3]) are written in a form similar to those 
for FLRW models, 



ds 2 = -dt 2 + A 2 



l + 



A, r T 



A 



dr 2 



1 — k(r)r' 



+ r 2 d£l\ 



(8) 



6) 

P(t,r) 



k(r) p (r) 



A 2 3A 3 
(Po^ 3 ),r 



3A 2 r 2 (Ar) t1 

The solution of Eqs. @ and (jSJ) is now written as 



afar) 
tfa r ) 



Qk(r) 

Pop) 
6k(r) 



1 — cos ^yj k(r) 7] 



7] 



1 



k(r) 



sin I W k(r) r\ 



+ U{r) 



(9) 
(10) 



(11) 
(12) 



where r] = rr = f t dt'/A(t',r) and A(t(r],r),r) = a(r),r). From the above definitions we can 
consider rj a generalized conformal time coordinate. 

Due to the freedom in the choice of the radial coordinate we can set the function po(r) to 
be a constant, po{r) = p n = constant, corresponding to the choice of coordinates in which 



M(r) oc r 3 , the so called |44| the FLRW gauge. We will adopt these coordinates in the rest of 
this paper. 



III. LOW REDSHIFT EXPANSION 

Since the geodesic equations are simpler if written in in terms of the variable t, it can be 
useful to find a perturbative solution of for R(t,r), expanding the Einstein's equations around 
the point (to, 0) corresponding to the central observer. 

We will expand the relevant functions in the following way: 



t(z) 
r(z) 
E(r) 
k(r) 



t + hz + t 2 z 2 + t 3 z 3 , 
r + r x z + r 2 z 2 + r 3 z 3 , 
&2r 2 + e 4 r 4 , 



k + k 2 r 2 , 



(13) 
(14) 
(15) 
(16) 
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where we have only even powers for E(r) and k(r) to ensure that the solution is analytical 
everywhere. In order to obtain a perturbative expansion for the solution it is convenient to 
introduce a new variable x and a function f(x) according to: 

M 

R = W\ j[x) ' (17) 

M p 

in terms of which the Einstein's equations corresponds to: 

\E\ , \ 

p=y (19) 

f(x) 2 -j-p = 0. (20) 



This choice of variables is different from 



44| , and in fact it leads to the simpler differential 



equation ( 120]) . Another difference is that the relations derived below in eq. (1241126]) for the 
coefficients of the expansion are valid for any sign of k(r) , while in 44| a two branches solution 
is used, requiring the necessity of the additional proof that the expansion is the same for each 
branch. 

The main advantage of re-writing the solution in this form is that we now have a constant 
coefficients differential equation in the single variable x, which is easier to Taylor expand. 
We can then solve this differential equation perturbatively around xq = x(to, 0) 

f(x) = /o + fi(x - x ) + ^(x - x ) 2 + ^(x - x ) 3 + ^(x - x ) 4 , (21) 



fi = \T+P / 2 = -4, (22) 
V Jo Jo 

Jo Jo 



f 2 r¥ f 2(3 /oP + 7) 

J3 - 73 U 75 • V Z6 I 



The above expansion of the solution has been obtained without any knowledge of the exact 
anaytical solution, but combining this latter with the definitions of f(x) we obtain: 



f(y) = -p(l-cosy), (24) 
x(y) = (-p) 3/2 (y - siny), y = Vkf], (25) 
fo = f(x ) = -po(l - cos J~k~ r]o), p = — = — (26) 



e 2 k Q 

which can be used to derive with an alternative method the coefficients of the expansion of 
f(x), and in particular to fix / . We can observe that fo, which encodes the initial conditions 



of the solution since all the other coefficients depend on it, depends only on the product 
\/^o?7o, implying a re-scaling symmetry on rj and k Q under which observational quantities are 
invariant. Applying recursively the derivative chain rule we can for example obtain the linear 
and second order terms using : 

df_ = d£dy = df_ f&A" 1 
dx dy dx dy \dy J 

d 2 f d ( df fdx 



(27) 



dx 2 dy \dy \dy 




(28) 
(29) 



It can be verified that the coefficients fi obtained solving perturbatively the differential equation 
fl20|) or using the above mentioned method are the same. These relations are also important to 
establish the a connection with the parameters appearing in the exact solution. 

In order to check the equivalence of the calculations done using the two different set of 
coordinates (r],r) and (t,r), it is also useful to define 



Vo 7f = sin (\/k r] 

V fc o 



+ t b (0) . (30) 



IV. GEODESICS EQUATIONS AND LUMINOSITY DISTANCE 

The geodesies equations are 1451 ] : 



dr jl + 2E(r(z)) 

dz (l + z)R'[r{z),t{z)Y { ' 

(M R'[r(z),t(r)] 
~dz~ ~ (l + z)R'[r(z),t(z)] ' 

(32) 

where the ' denotes the derivative respect to r and the dot ' the derivative respect to t. These 
equations are derived from the definition of redshift and by following the evolution of a short 
time interval along the null geodesic T(r). The r.h.s. can be evaluated from the perturbative 
solution for R(t, r) in terms of f(x). Expanding in powers of z we can then reduce the solution 
of this system of partial differential equations to solving a system of linear algebraic equations, 
where the unknowns are the expansion coefficients: 

G(z) = — = G + G lZ + G 2 z 2 , (33) 



dt 

L{z) = — = Lo + L lZ + L 2 z 2 , (34) 
az 

r(z) = r + T\z + r 2 z 2 + r 3 z 3 , (35) 

t(z) =t + t lZ + t 2 z 2 , (36) 

n = G , 2r 2 = d, 3r 3 = G 2 , ^ = L , 2t 2 = L x . (37) 

We will assume also assume a homogeneous big-bang, = In the above equations we are 
expanding all the quantities to the order necessary to expand the luminosity distance to the 
third order. For the geodesic equations we get: 



/ \ / fo v / 7o(/oPo + 3) 2 , . 

r[Z) ~ V 2e 2Po (/oPo + 2) Z 2v/2^(/ R, + ^ K } 



1 
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24eip (/oPo + 2) 3 
3v / 2e 4 / 3/2 v ^(/ Po + 2) 3/2 



e\ (2po v / 2/oe 2 Po(/oPo + 2) (3/ 2 po + 12/ + 15po) + 216e 4 to(/oPo + 2] 



t(z) = t - f° /2 P°y/W° z + /0 /2 P0V / ^( 2 /0P0 + 5 ) ^2 f39l 

W ° 12>/2e§v7S^+2 24y2e|(/ p + 2) 3 / 2 ' 1 ' 

Using eqs.( )24ll26|30|) it is more elegant and convenient to re-express our results in terms of 
r)o ,ki, and po which gives the following formula for the luminosity distance Dl(z) : 



D L (z) = (1 + z) 2 R(t(z), r(z)) = (1 + z) 2 r(z) 2 a( V (z),r(z)), (40) 

D L (z) = D\z + D\z 2 + D 3 L z 3 + ... 

B 3 p 
fco /2 (3£ 2 + 3)' 
B 3 (B 2 -l)p 
12kl /2 (B 2 + 1) ' 

B 2 p (k 2 (B 4 - 1) + 8k 2 B 2 X + 2k 2 (3B 4 X - 9B 3 + 8B 2 X - 9B + 9X)) 



£)3 = 

L 2Akl /2 (B 2 + 1) 

X =^rj /2; 5 = tan(X). (41) 

In the rest of the paper we will not give the formulae in terms of r/o and trigonometric func- 
tions, since they are rather complicated. So far we have in fact worked in terms of functions and 
parameters which depend on the coordinate choice, but since they are not directly observable 
quantities, it is more convenient to introduce : 



B 2 Po /, Q v 

ao = 3k (B 2 + iy (42) 
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Hlt = !M^!±il, (43 ) 



\ (B 2 + l) , (44) 

where we have used the following definitions 



go 2 



a = a(r) , 0) , (45) 

H ° = < 46 > 

a(t o ,0) 

afo, 0)d(to, 0) 
a(t ,0) 2 

From the equations for go we can see that it must be always positive. As we will see 
later, under the assumption of central smoothness we have adopted we have that g = go PP , in 
agreement with previous studies ^6], which proved that the apparent deceleration parameter 
cannot be negative for a smooth LTB solution. 

The formula for the luminosity distance takes now the form: 

D\ = ^, (48) 
-"0 



D 2 - 

" 2H n ■ 



D s = (go - ljgta + 3K2% 



, ^ arctan(V2q -l) 

2 Wo + 1) — 



L ~ 2H 2H (2q -iy 
where we have introduced the dimensionless parameter K 2 = k 2 / '(aoH ) 4 . The inhomogeneities 
effects show only from the third order coefficient D\ because we have not included odd powers 
in the expansion of k(r) in order to satisfy smoothness conditions, and in the {K 2 = 0, go — 1/2} 
limit we recover the case of a flat matter dominated FLRW Universe, i.e. Qm — 1- 

From the definitions of H , g and a we get the following relations with the parameters 
appearing in the exact solution: 



X = arctan y / 2g - 1, (49) 
k = 4ff 2 (2g -i), (50) 
p = Qa 3 H 2 q . (51) 

These formulas are valid for both negative and positive ko by analytical continuation, and 
provide a very compact and insightful picture of the physics of the problem. Eq. (|50p gives 
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in fact the direct relation between q and the value of ko, and implies that q > 1/2 only for 
positive ko. 



V. APPARENT COSMOLOGICAL OBSERVABLES 



The following relations apply to a flat FLRW model 



H app (z) 



dz \ 1 + z 



(52) 



Q app {z) = ± (y^j) = (H^(z))' 1 , (53) 

Apparent observables are then defined as those obtained by applying the same relations above 
to the luminosity distance D^ TB (z) obtained for a central observer in a LTB model. Another 
investigation of the effects of inhomogeneities on apparent cosmological observables in presence 
of a cosmological constant can be found in j^], where is shown that the inhomogenity could 
introduce the illusion of a red-shift dependent equation of state, instead of a cosmological 
constant. It can be easy verified that 

q app (0) = qo PP = 9b, (55) 
H app (0) = H% pp = H 0} (56) 

where it is important to observe that q and qQ PP are defined independently in eq. (H7|) and 

ity between q and qQ PP only holds under the assumption that 



eq. (j54]l respectively. The equa 
k'(0) = t' b (0) = as shown in j48j. This can be considered a natural consequence of the fact 
that for a smooth LTB model the apparent deceleration parameter q^ at the center should 
coincide, by continuity, with go- As expected, in the centrally smooth case we are considering 
qQ PP is always positive because of eq. (I44|) . This is an alternative j^l and easier prove of the 
fact that a smooth LTB space has a positive central apparent deceleration parameter: 

qo PP > 0. (57) 

After substituting the expansion of the luminosity distance in the definition of apparent cos- 
mological deceleration given in eq. (l54"j) we get 



\ (B 2 + l) = g , (59) 
apP k 2 {B 2 -1)B 3 + 9k 2 (B 2 + 1) (B 2 X -3B + 3X) 

1l = 7^TT> • ( 60 ) 



As it can be seen, since we consider a centrally smooth model, the effects of the inhomogeneity 
show only in the linear order coefficient q1 PP , but this will be enough to derive the general leading 
order conditions for a negative deceleration at a red-shift different from zero. 



VI. CONDITIONS FOR A NEGATIVE APPARENT COSMOLOGICAL DECELER- 
ATION 

From the expression derived in the previous section we can deduce the leading order condition 
to have a negative q app (z n ) < at a given redshift z n 

q app < -~ql pp - (61) 
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FIG. 1: q\ PP {qo) is plotted for different values K 2 = ^/{oqHq) and for < qo < 1. The small dashed 
line corresponds to K 2 = 10, the medium dashed line corresponds to K 2 = and the long dashed line 
corresponds to K 2 = — 10. 



Since gg PP is always positive, we will focus here on the linear order correction, and after intro- 
ducing the dimensionless parameter K 2 = ^/((IqHq) and using eq. (15o]l . we can express the first 
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FIG. 2: Qi PP (qo) is plotted for different positive values of K% = ^/C^q-Hq) = 1)2,4,8 and for 
< (fci < 1. The straight lines correspond to q1 PP = —\% for different values of z, i.e. with different 
slopes. Negative values of q app (z n ) correspond, from eq. (f6Tj) . to the area of the (qo, q\ PP ) plane below 
the corresponding straight line with equation q°^ p = ——qo- This implies that, keeping K2 fixed, 
for smaller redshifts, i.e lower slope lines, the maximum value of qo which gives a negative q app (z) 
decreases. We can also see that at the some redshift, i.e. for the same straight line, higher values 
of K2 give a larger maximum value of qo necessary to satisfy the condition for a negative q app (z). 
Increasingly Higher values of K2 correspond to curves with increasingly lower values of qi PP (qo)- 

order coefficient in the following form: 

From this last formula we can see that the only physically relevant parameters are the dimen- 
sionless quantities go an d K 2 . From the figs. 1-2 representing q\ pp {qo) we can see that larger 
K2 allow q app (z) to be negative for smaller maximum values of go- Because of eq. (!50l) go and 
K are in fact proportional, and so a smaller value of K is compensated by a faster growth of 
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k(r), which is the same reason why negative K 2 , corresponding to a decreasing k(r), requires 
larger values of go- The straight lines drawn in the fig. [2] have slope — -, and from eq. flBTj) a 
negative q app (z) correspond to the area of the (go, Qi PP ) plane below those critical lines. This 
implies that, keeping K 2 fixed, for smaller redshift, i.e lower slope lines, the maximum value of 
go which gives a negative q app (z) decreases. From the same figure we can see that at the same 
redshift, i.e. for the same straight line, higher values of K 2 give a larger maximum value of go 
necessary to satisfy the condition for a negative q app (z). 

VII. COORDINATE INDEPENDENT CONDITIONS 

So far we have derived the conditions for a negative apparent deceleration parameter under 

3 

the choice of the radial coordinate in which M(r) = but it would be more useful to derive 
a set of coordinate independent conditions. This can be achieved by considering the expansion 
of the relevant functions in terms of the red-shift: 

k(z) = k Q + k{z + k z 2 z 2 + .. (63) 

After substituting the expansion for r(z) in k(r) we get 

k(z) = k(r(z)) = k + (fcin)z + ^(2k 2 rl + hr 2 )z 2 , (64) 



from which we obtain the relations 



fci = — , (65) 

n 



— 2Hn + k z r 2 

k 2 = — 2 2 ; f 1 > (66) 



which imposing the smoothness condition reduce to 

h = 0, (67) 

k 2 = M . (68) 
'1 

In the above expression the coefficients ki,ri are the same as the ones defined in equations 

We can now substitute the above relations in the formulae obtained in the previous section, 
finally getting an expression for q app (z) which only depends on the coordinate independent 
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coefficient of the redshift expansion K\: 

<Ko(K* - Kl + 27 Kl + 27 K K*) + 



q°*P( z ) = q * pp + q f* z + .. = + ' 



2 2KT 



-9(3 + AK + K$)Kl arctan JK 



z + .., (69) 



where we have expressed the final result in terms of the dimensionless parameters 

K = k (a H )- 2 = (2q - 1) , (70) 
K z 2 = k z 2 (a H y 2 . (71) 

The condition for a positive acceleration remains the same as in eq.f lSTTh 

It should be noted that while the expansion of r(z) is different for different choices of the 
radial coordinate, that for for k(z) is coordinate independent, implying that the above expres- 
sion for q app (z) can be applied to any LTB model. It can be checked explicitly in fact that 
for example adopting the light cone gauge the expansion for r(z) would be different, but the 
formula ( |69l) would be the same since it is a physical observable. 
The redshift expansion of k(z) would also be the same since 

k fr [r fr (z)} = k fr [r fr (r lc (z)} = k lc [r lc (z)} = k(z) , (72) 
k lc [r lc ] = k fr [r fr (r lc )), (73) 

where for clarity we are denoting with r; c the light cone gauge radial coordinate and rj r the 
FLRW gauge radial coordinate. A direct check would require to derive explicitly a local expan- 
sion of the coordinate transformation between r/ r (r; c ) which can be obtained by imposing the 
condition: 

M' c K(r>)] = M^[r /r ] = ^±, (74) 
ric(r fr ) = r[ c r fr + (r>) 2 + r l 3 c (r fr ) 3 + .. (75) 

This is a rather cumbersome procedure, so we will not report it here, but simply observe that 
our final result in eq.( )69l is coordinate independent because it involves the coefficients of the 
expansion of k(z) and not of k(r). A graphical representation in the plane Ko,qi PP of the 
condition for a positive acceleration is give in Fig. 3. 



13 



app 
?1 
0.4 r 




FIG. 3: ql pp is plotted for different values of X| = {-0.2, -0.1,0.1,0.2} and for -1 < K < 1. The 
small dashed line corresponds to z = 0.01, the medium dashed line corresponds to z = 0.1 and the long 
dashed line corresponds to z = 5. Negative values of q app (z n ) correspond, from eq.(|61|) . to the area of 
the (Ko,q1 PP ) plane below the corresponding straight lines with equation q^ pp = —j-Qo = — — 1+ 2 ^° • 
Higher values of K% correspond to lower values of q± . 

VIII. CONCLUSIONS 

We have derived the low-redshift expansion for the apparent cosmological deceleration pa- 
rameter q app (z) for a central observer at the center of a centrally smooth spherically symmetric 
matter inhomogeneity. The same results have been derived with two different methods, one 
based on the existence of an exact solution in terms of coordinates {f],r), and the other using 
a local perturbative expansion of the solution around the center using cosmic time. We have 
then have applied it to study the conditions for a negative value of the apparent deceleration 
parameter q app (z) at non zero redshift, showing how the maximum allowed value of go decreases 
as the redshift decreases, and it decreases as K 2 decreases. We have finally derived a coordi- 
nate independent expression for the apparent deceleration parameter, based on the red-shift 
expansion of the relevant functions. 

The method we used to expand the solution of the Einstein's equation has the advantage 
of using cosmic time and not the parametric analytical solution, allowing to use it for other 
applications which require space-like foliation of space-time such as volume averaging. 

In the future it will be interesting to apply our analytical methods to study other relevant 
cosmological observables for a central observers in LTB spaces, in order to better understand 
the local effects effects of large scale inhomogeneities. 
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Appendix A: Alternative method 



Here we adopt the same method developed in 49| to find the null geodesic equation in the 
coordinates {rj,t), to find a local expansion of the solution around z = corresponding to 
(t , 0) = (i] ,t), where t = t(f]Q,r). The luminosity distance for a central observer in a LTB 
space can be written as 

D L (z) = (1 + zfR (t(z), r(z)) = (1 + zfr{z)a ( V (z), r(z)) , (1) 

where (t(z),r(z)^j or {[r]{z), r(z)j is the solution of the radial geodesic equation as a function 
of the redshift in different coordinates. 
The radial null geodesic equation is 

dT{r) -R r (t,r) 

— — = / T r ,r ; f(t,r) = —= 2 
dr <'l + 2E(r) 



where T(r) is the time coordinate along the null radial geodesic as a function of the the coor- 
dinate r. From the analytical solution, we can write 

T(r) = t(U(r),r), (3) 
dT{r) _ dt dU(r) + dt 
dr dr] dr dr ' 

where U(r) is the r] coordinate along the null geodesic as a function of the the radial coordinate 

r. It is then possible to write the geodesic equations for the coordinates (r/,r), 

drj d r t(rj,r) — F(r],r) 



dz (1 + z)d v F(r],r) 
dr a {Vi r ) 



dz (1 + z)d rj F(r), r) 



p(f],r), (5) 
q(r],r), (6) 



F(r/,r) = [d r (a(r],r)r) + d v (a(T],r)r)d r T]} , (7) 

Jl — k(r)r 2 
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where rj = U(r(z)) and F(r),r) = f(t(j],r),r). It is important to observe that the functions 
p, q, F have an explicit analytical form which can be obtained from a(r], r) and t(r), r). 

In order to take advantage of the fully analytical expression of the equation of the radial null 
geodesies we can follow the same procedure we adopted for the coordinate (t,r), and expand 
the geodesic equations solution as: 



T](z) = T/o + TJxZ + 7] 2 Z + ... 



(9) 



The result for rj(z) is 



2 arctan ( \/2q — 1 1 



V2qo - 1 (Ho (4<? - 3<?o + 1) - 3A; 2 ) + 2k 2 (q + 1) arctan U2q - 1 



z 2 



2H$(2q - 1)^ 

while r(z) is the same as the one obtained in the previous sections since the coordinate r does 
not change. 

It can be checked that the two methods to calculate the geodesies are equivalent by substi- 
tuting in the exact solution t(r),r) 

t(z)=t(r]{z),r{z)), (11) 

the expressions r](z) and r(z) and comparing with the results for t(z) from the previous sections. 
Such a substitution gives exactly the some expression for t(z), proving that the two methods 
to obtain the red-shift expansion of the solution of the geodesies equations are equivalent. 
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